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基于 谱 对 称 的 形状 配 准 方法 
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摘 要 : 针对 现 有 三 维 形状 配 准 方法 中 存在 左右 翻转 的 错误 匹配 问题 ， 提 出 了 基于 内 蕴 对 称 特 征 检测 的 高 效 形状 配 准 
算法 。 首 先 ， 通 过 热 核 与 几何 约束 ， 构 建 模型 的 内 蕴 自 对 称 点 对 ; 其 次 ， 基 于 谱 误 入 特征 空间 分 析 ， 提 取 模 型 的 内 区 


对 称 平面 ， 并 依据 模型 表面 法 向 量 ， 有 效 识 别 模型 的 左右 结构 属性 ; 然后 ， 根 据 内 列 对 称 点 对 ， 获 得 模型 的 一 致 性 谱 
对 称 结构 描述 ; 最 后 ， 引 入 一 致 性 点 漂移 算法 (CPD)， 实 现 基于 谱 对 称 的 非 刚 性 模型 的 形状 配 准 ， 有 效 避 免 了 模型 配 准 
中 的 左右 结构 翻转 问题 。 实 验 进一步 论证 了 这 种 方法 不 仅 有 效 提 高 模型 匹配 的 效率 ， 而 且 能 有 效 识 别 同类 模型 的 结构 


特征 ， 对 于 非 刚 性 模型 的 配 准 有 具有 较 强 的 鲁 棒 性 。 
关键 词 : 形状 配 准 ; 对 称 性 检测 ; 谱 映 射 ; 一 致 性 点 漂移 算法 
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Shape registration method based on spectral symmetry 


Han Li, Li Dan, Liu Shuning, Chen Siyu, Pu Jingyu, Tangdii 
(School of Computer & Information Technology Liaoning Normal University, Dalian Liaoning 116081, China) 


Abstract: Aiming at the issue of symmetric flips in the process of 3D shape registration, this paper developed an efficient shape 
registration algorithm based on intrinsic symmetric feature detection. Firstly, it constructed intrinsic symmetric point pairs of 
the model by heat kernel signature (HKS) and geometric constraints. Secondly, based on the spectral embedding space analysis, 
it extracted the intrinsic symmetric plane of the model and effectively identified the symmetrical properties of the model 
according to the model surface normal vector, getting intrinsic symmetry point pair. Therefore it presented the consistent spectral 
symmetry structure of the model. Finally, combining the coherent point drift method, it implemented the shape registration of 


non-rigid model based on spectral symmetry. The experimental results show that the matching method is efficient and robust to 


the non-rigid deformable shape matching. Moreover, the structural features in same category models also are effectively 


— identified. 
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Kazhdan 等 人 所 研究 了 一 种 对 称 变换 ， 用 于 度量 二 维 形状 关 了 
平面 内 任意 直线 反射 对 称 的 可 能 性 。 这 种 方法 后 来 被 Podolak 
在 几何 模型 的 形状 分 析 中 ， 借 助 对 称 性 可 以 消除 模型 中 的 等 人 BJ 扩展 应 用 于 三 维 几 何 模型 ， 提 出 了 平面 反射 对 称 变换 
元 余 信息 ， 可 以 完成 对 模型 的 高 层次 的 表示 和 分 析 ， 亦 可 以 完 (planar reflective symmetry transform, PRST)» Rustamov!“ 
成 对 模型 的 对 称 化 、 残 损 模 型 修补 等 模型 优化 处 理 。 因 此 ， 研 体 非 对 称 性 的 空间 分 布 结合 PRST 中 ， 得 到 了 增强 的 PRST 
究 三 维 模型 的 对 称 性 具有 重要 的 理论 与 工程 意义 。 用 于 高 效 的 形状 检索 。Martine 等 人 中 通过 分 析 球 面 调和 系数 定 
按照 对 称 性 定义 所 依赖 的 几何 属性 可 分 为 外 蕴 对 称 和 内 蕴 义 了 一 种 广义 矩 ， 并 基于 此 提出 一 种 恢复 三 维 形 状 全 局 对 称 性 
对 称 ; 按照 对 称 在 三 维 模型 表面 的 分 布 ， 又 可 以 分 为 全 局 对 称 的 方法 。 在 局 部 外 列 对 称 研究 中 ,Thrun 等 人 中 提出 了 一 种 分 级 
和 局 部 对 称 。 的 产生 一 测试 过 程 用 于 局 部 对 称 检测 。Simari 等 人 中 研究 了 局 
大 多 数 对 称 性 分 析 工 作 集中 于 全 局 外 列 对 称 的 检测 和 分 类 。 ”部 近似 平面 反射 对 称 的 检测 问题 ， 并 为 三 维 形状 定义 了 一 种 基 
Zabrodsky 等 人 [定义 了 一 个 变换 的 对 称 性 度量 , 用 于 衡量 一 个 于 平面 对 称 子 部 分 的 层次 表示 。Mitra 等 人 四 提出 了 一 种 变换 空 
形状 到 它 的 最 近 对 称 形状 的 距离 ， 实 现 了 近似 对 称 检测 。 间 选 举 (yoting) 与 聚 类 (clustering) 的 方法 来 检测 局 部 近似 对 称 ， 
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并 进 


步 研究 了 三 


其 外 列 对 称 四 。 


在 内 区 对 称 检 测 方 


fl, H 


Lipman $ ATONE H 
检测 方法 ， 该 方法 为 输入 点 云 构 建 对 称 对 应 和 矩阵， 然后 利用 
上 定义 了 与 扩散 距离 (diffusion 


矩阵 的 频谱 特性 ， 
distance) 密 切 相 关 的 对 称 分 
Riklin-Raviv 等 人 03 针 对 流 形 曲 面 上 的 一 个 采样 点 集 ， 提 出 了 
自 底 向 上 搜索 该 点 集 与 其 自身 的 等 距 变 换 。 全 
point signatures, GPS) 可 以 将 形状 的 全 局 内 区 对 称 变换 为 高 维 空 


佳 模型 的 对 称 化 ， 即 对 一 个 模型 施加 形变 使 


前 绝 大 多 数 方法 仅 考 虑 全 局 内 蕴 对 称 。 


BB 了 一 种 三 维 点 云 
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维 曲 面 


的 对 称 群 轨道 (orbits) 的 


方程 描 ; 
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_ Ou(x,t) 
ôt 
Hp: Ay XM EH Laplace- Beltrami 算 子 ， 如 果 M 有 边 


Ayux,t) = (1) 
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分 解 距离 (symmetry factored distance) . 


界 , 使 U WE Dirichlet 边界 条 件 , 即 对 于 流 形 M 上 所 有 点 并 
MAE t, Hpi u(x.) = 0 ; u(x,t) 为 在 时 间 t 时 刻 位 置 
上 的 温度 。 给 定 一 个 初始 热源 f  ， 热 传导 方程 在 时 间 t+ 上 的 解 


局 点 签名 (global 


间 上 的 全 局 外 蕴 对 称 na44。 类似 地 , Chertok 等 人 U5 研究 了 基于 
频谱 分 析 的 对 称 检测 方法 。Kim 等 人 09 提 出 了 基于 Mobius 变 
换 的 全 局 内 蕴 对 称 检测 方法 。Zhang 等 人 [7 提出 了 一 种 形状 指 
述 算法 (symmetry robust descriptor，SRD)， 可 以 鲁 棒 性 的 区 分 两 
对 有 角度 标志 的 点 ， 甚 至 是 对 称 或 者 近似 对 称 的 点 。Zheng 等 
T 和 人 09 提出 对 形状 的 内 莹 骨架 的 对 称 来 进行 研究 ， 增 强 形状 的 自 
= 我 对 称 性 。 
CA 近年 来 , 非 刚性 形状 对 应 研究 取得 了 显著 的 进展 Ht9221。 其 中 
N 最 具有 代表 性 的 方法 有 IsomapP23 和 Laplacian 特征 映射 29， 其 
将 一 个 内 在 的 形状 匹配 问题 转化 为 一 个 高 维特 征 空间 中 的 外 在 
ef 距离 度量 问题 .但 是 由 于 高 维特 征 空间 的 符号 翻转 以 及 重复 性 ， 
© 导致 在 对 称 检测 及 对 应 匹配 中 出 现 左右 结构 翻转 问题 。 
= 文献 [25] 提 出 对 应 对 称 性 感知 柑 入 ， 与 以 前 将 表面 嵌入 高 
O 维 空间 的 其 入 方法 不 同 ， 而 是 基于 低 维 嵌入 的 。 该 方法 通过 在 
N 两 个 没有 反射 的 嵌入 表面 之 间 找 到 3D 刚体 变换 来 解决 左右 翻 
DO 转 问 题 ， 但 是 该 方法 局 限于 耗 时 的 迭代 求解 以 及 不 能 处 理 旋转 
DC 对 称 和 重复 结构 。 
= 本 文 基于 谱 图 理论 ， 提 出 了 内 蕴 自 对 称 平面 的 检测 方法 ， 
= 有 效 揭示 模型 的 左右 及 形状 结构 ， 并 引入 谱 对 称 的 概念 ， 实 现 
O 高效 精准 的 形状 匹配 。 首 先 ， 通 过 热 核 与 几何 约束 分 析 模型 的 


自 对 称 属性 ， 提 


HDL 
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自 对 称 点 对 ; 其 次 ， 
型 内 莹 对 称 平面 的 计算 ， 构 建 模 型 的 结构 特 怕 


对 称 点 对 提取 模型 的 一 致 性 谱 对 称 稀 下 C 


方 法 ， 实现 高 效 精准 
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BAR ACH 


本 文 方法 不 仅 有 效 避 免 J 
j 且 稀疏 化 和 结构 化 的 谱 对 


翻转 问题 ，T 
模型 的 配 准 效率 。 


模型 配 准 


与 分 析 。 
中 出 现 的 左右 以 及 前 后 


基于 低 维 谱 典 入 进行 模 


Ea, FET A Bi 
HB; 最 后 ,引入 CPD”? 


可 以 通过 热 算 子 H, 得 到 ， 即 
u(x,t) =H, f (2) 
对 于 任意 流行 M FERS K, x,y): R xMxM OR 
使 得 


u(x,t) = | K(x,y)f Oy)dy 


M 


(3) 


其 中 : 


d, 是 yeM 的 体积 形式 ; 满足 该 方程 的 最 小 值 


Ki(x,y) 称 为 热 核 ， 可 以 看 做 经 过 时 间 t 后 点 xX 处 单位 热量 传 


递 到 点 》 的 热量 总 和 。 热 核 具 有 以 下 的 热 核 分 解 : 


K,(x,y) = 》e OO) (4) 


i=0 
其 中 :元 AAMA Laplace — Belrtami 算 子 的 第 i 个 特征 
值 和 对 应 的 特征 向 量 。 
根据 Coifman 提出 的 扩散 距离 框架 ， 可 度量 
两 数据 点 之 间 的 相似 度 。 对 模型 M 
距离 : 


标 空间 中 
， 给 定时 间 z 可 定义 扩散 


D, (x,y? = P z)-K,(y, 2) dz (5) 


K) ARAB, 7 PULA EB KA H 
数 点 集 之 间 的 扩散 距离 等 价 于 目标 空间 的 欧 氏 距离 。 
热 核 特征 记录 经 过 一 段 时 间 后 点 的 热量 扩散 到 其 他 
点 ， 可 以 得 到 模型 上 该 点 邻 域内 的 几何 信息 。 显 然 ， 热 核 特征 
值 由 时 间 域 所 主导 ， 有 着 时 间 域 上 的 多 尺度 特性 。 某 一 点 X 
的 热 核 特征 值 可 看 做 时 间 z 变化 的 时 间 序 列 : 
K,O, xi) = (Kx, xX;), Kyo (Xi, Xi),..., Ky (Xi, Xi)) 

1 所 示 为 上 = 45 时 刻 不 同 姿态 模型 上 的 热 核 分 布 情 
况 。 图 中 由 红色 到 蓝 色 ( 见 电子 版 〉 表 示 热 核 特征 值 由 大 到 小 


标 空 间 ， 


称 表示 更 进一步 提高 J 


1 ”基于 热 核 与 几何 约束 的 内 蕴 自 对 称 分 析 


本 文 利用 热 核 特征 的 内 蕴 等 距 不 变性 ， 进 一 步 结合 几何 形状 约 
束 ， 提 出 了 高 效 的 内 药 自 对 称 分 析 方 法 。 

1.1 热 核 特征 

热 核 描述 符 (HKS)P 9 是 一 种 基于 热 扩散 属性 描述 模型 表面 几 
何 特征 的 描述 符 , 在 保持 模型 内 蕴 几 何 属性 的 同时 ,构造 简洁 ， 

对 噪声 不 敏感 。 

MT ABR Site M (可 能 存在 边界 ), 其 热 扩 散 过 程 可 以 用 热 


x 


变化 ， 可 见 具 有 相似 对 应 结构 关系 的 子 部 分 其 热 核 函 数值 相 


Tk 


3D 模型 的 热 核 分 布 


图 1 


然而 热 核 函数 却 不 能 精确 识别 模型 的 形状 细节 ， 例 如 人 的 模型 
中 ， 手 和 脚 的 热 核 值 极为 相似 ， 若 仅仅 通过 热 核 特征 去 分 析 模 
型 的 对 称 属性 ， 就 会 出 现 错误 的 手脚 对 称 特征 。 因 此 ， 本 文 进 


录用 稿 


一 步 融 合 模 型 的 几何 约束 ， 突 出 模型 形状 细节 ， 实 现 稳 

确 的 自 对 称 形状 检测 。 

1.2 融合 热 核 与 几何 约束 的 内 蕴 自 对 称 提取 

测 地 线 距 离 可 获得 高 维 观 测 空 间 与 低 维 流 形 结构 在 局 部 意义 下 
的 内 在 几何 分 布 结构 ， 能 够 有 效 地 揭示 流 形 空间 的 内 在 特征 。 
本 文 融合 热 核 的 等 距 不 变性 以 及 测 地 线 距 离 的 内 萄 结构 性 ， 提 
出 了 健壮 的 三 维 自 对 称 形 状 分 析 算 法 。 
本 文采 用 无 向 带 权 图 CV, E) 来 表示 模型 MM ， 


其 中 ， v EV 


X 


为 模型 表面 项 点 ，E 代表 每 一 对 顶点 Vi 和 Vj 的 连接 边 ， 则 


G 的 邻接 矩阵 ACG) = (aj) EA 


(6) 


Hh: woy) 为 对 应 边 Cy 的 权 值 。 为 了 突出 模型 顶点 间 的 局 


H 邻 域 关 系 ， X) 


高 斯 核 函数 (Gaussian kernel) 来 构建 图 G 的 相 


WIRE REE Ws = (Wij) ，W。e R” 每 对 顶点 的 相似 程度 用 点 之 


间 边 的 权 值 wj s[0,]] 表示 。 


-dist(v; Vj y 


Jep: distav) 表示 顶点 Vi; 与 Vj 之 间 的 测 地 线 距离 ，0 为 


高 斯 分 布 的 尺度 参数 。 
首先 , 计算 模型 M 的 热 核 特 征 , 点 的 热 核 特 征 向 量 构 成 了 一 个 


规范 正 交 系 ， 模 型 M 可 用 该 集合 的 基底 {9 ),,| 表示 : 


Kj = > ago; K, 一 4,9; (8) 


i>1 j2l 


n 为 模型 M 的 顶点 个 数 ; 为 


Jp: Ka Ky k 维 向 量 ; 


所 选取 的 特征 值 的 个 数 。 
其 次 ， 通 过 度量 热 核 特征 判断 模型 的 自 对 称 属 性 ， 根 据 式 (4) 
和 “5)， 顶 点 之 间 的 热 扩 散 距 离 定义 为 


di tp Xi Xj) 一 


Rt) Xt) 


d i Ye” 


) dlogt)’* (9) 
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Fe: KC) 为 热 核 ，> ,e “为 1 时 刻 热 核 的 迹 。 可 以 描述 


数据 点 之 间 的 局 部 连通 性 ， 在 对 数 上 均匀 采样 按 比例 缩放 ， 得 


到 一 个 向 量 来 表示 HKS 每 个 点 。 表 1 为 特征 向 量 Ks; M Kaj 


构成 hnXn 热 核 特 征 相 似 和 矩阵 。 
表 1 TURE AT REAL EAE BE 
K, Ka K,, K, 
K, d(1,1) d(1,2) d(1,3) d(l,n) 
Ky, d(2,1) d(2,2) d(2,3) d(2,n) 
K,, d(3,1) d(3,2) d(3,3) d(3,n) 
K,, d(n,1) d(n,2) d(n,3) d(n,n) 


JE PE HE HI BE — ITAR = HT TOL AA Sp SH ft T AA S FAK A RE 


( 亦 为 对 称 度 )。 


进而 基于 获得 的 热 核 对 称 关系 矩阵 dgxs Vio Vj) 与 测 地 线 距离 


相似 度 和 矩阵 W。 ， 构 建 几何 约束 的 自 对 称 检测 矩阵 : 


Dym (Vi, Vj) = Ad pgs Vi Vj) + PW, (vi,v)) 


其 中 : a+ B=1 .实验 中 ,为 了 突出 模型 的 局 


部 几何 属性 特征 ， 


设置 权 值 w =0.87 ， B= 0.13 。 Domi Vj )) 表示 模型 上 的 顶点 


vi 与 Vj 之 间 的 对 称 程度 , 遍历 整个 矩阵 找 出 最 小 值 ,获得 最 优 


的 对 称 点 对 : 


min(D pym (V ni? Vaj)» Il ii) (10) 


其 中 : 7 为 模型 点 的 个 数 ，TT; Ati Ae, wert eae 直 


到 完全 遍历 整个 矩阵 。 其 算法 如 下 : 
算法 1 基于 热 核 特征 与 几何 约束 的 自 对 称 检测 算法 
输入 : ”三 角 网 格 模型 M 。 


输出 : 模型 的 对 称 点 对 集 y 。 


了 计算 三 角 网 格 模型 的 Laplace-beltrami 算 子 ; 
b) 计 算 Laplace-beltrami 算 子 的 特征 值 4 和 特征 向 量 9; ; 


9 计算 时 间 te [inin,tnax] 范围 内 的 模型 热 核 特征 Ki(.,.) ; 


中 计算 顶点 间 的 扩散 距离 djyxs 以 及 测 地 线 距离 W。 ,构建 


自 对 称 检测 矩阵 Dam ’ Dyn = ad ygs T Bw, A 


e 通 过 遍历 对 称 对 应 矩阵 Dsym 获取 最 终 对 称 点 对 集合 
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T4t wt 


对 于 不 同 姿 态 的 几何 模 


2 


用 的 
关 重 
谱 对 
2.1 


图 2 不 同 模型 的 


图 2 为 基于 热 核 与 几何 约束 的 模型 自 


基于 谱 对 称 的 模型 配 准 算法 


三 维 


模型 的 配 准 


U pear 


对 称 检测 


对 称 检测 结 H 


与 形状 分 析 是 模型 检索 、 压 缩 、 重 建 等 应 


核心 内 容 ， 高 效 、 高 精 


要 的 作用 。 本 文 在 内 药 


的 形状 匹配 方法 研究 更 是 发 挥 着 至 
对 称 检测 的 基础 上 ， 提 出 了 基于 


称 的 非 刚性 模型 形状 配 准 分 析 方法 。 
2D RRA 
基于 Laplacian 算 子 的 谱 


央 射 将 


有 效 


wee. 


征 ， 


D 谱 特 征 


高 维 数据 


HRA BIR AE Ae 


间 ’ 


保持 模型 的 内 区 形状 特性 ， 被 广泛 应 用 于 模型 的 分 类 与 匹 


30。 文 献 [13] 指 出 越 大 的 特征 向 量 


包含 越 为 显著 的 形状 特 


因此 本 文采 用 拉 普 拉 
描述 ， 其 不 仅 直观 再 


现 模型 


免 高 


维特 征 空间 的 不 对 称 性 及 扭曲 性 。 


斯 映射 中 最 大 及 次 大 特征 


向 量 ， 构造 


4 的 形状 结构 ,而且 有 效 避 


但 是 拉 普 拉 


斯 特征 


映射 


所 提 
的 测 


供 的 嵌入 在 一 些 意识 中 只 保持 了 局 
试点 描述 不 充分 。 


F) 


是 | 


特殊 的 


斯 保 


对 于 
发 生 
均 发 
必然 


入 内 
称 


持 局 部 信息 的 最 佳 性 。 
如 图 3 所 示 ， 同 一 类 模型 的 二 维 谱 
姿态 与 局 部 变化 具有 较 强 的 稳定 性 。 


部 信息 的 最 


竺 性， 对 新 


对 称 感知 来 进行 形状 匹配 ， 在 进行 模型 的 配 准时 也 
局 部 配 准 (对 称 谱 域 配 准 )， 所 以 充分 利用 了 


拉 普 拉 


具有 结构 相似 性 ， 
而 谱 特 征 空间 极 易 


符号 偏转 现象 ， 
生 了 结构 翻转 ， 若 基于 谱 特 征 
造成 模型 匹配 过 程 中 的 左右 混淆 错误 。 
本 文 结 合 二 维 谱 图 上 
区 对称 平 面 ， 有 效 


图 3 中 不 同 姿态 的 同类 模型 在 二 维 谱 
描述 来 识别 模型 的 对 应 性 ， 


有 保持 模型 内 蕴 形 ; 
判断 模型 的 左右 结构 ， 


和 述 ， 实 现 基于 谱 对 称 的 模型 配 准 与 相似 怕 


域 中 


大 特性 的 优势 ， 引 
建立 优化 的 谱 对 
形状 分 析 。 


Jee 


CA 


图 3 不 同 姿态 及 不 同 种 类 模型 的 


维 谱 映射 


2.2 AAW 
本 文 依据 热 核 与 几何 约束 提取 的 


在 谱 域 中 获 


Q= (Qj0, Qjie 


I 
al, O= > 
g +Q 


面 拟 合 


(a) 谱 域 中 的 对 称 中 点 ”(b) 空 域 中 的 对 称 


They 


ARE A TARF 
结构 的 表示 方法 。 
设 模 型 的 内 莹 对 称 平面 P 


模型 


距离 d,: 


Di 


平面 进行 模型 检测 ， 


得 对 称 点 对 的 特征 


图 4-a 的 对 应 红 


平面 


Ree Myy); 


向 量 描述 工 = 


对 称 点 


Cio Ti. oTi) , 


Ox) ; 进而 衡量 特 


-0| 


依据 最 小 二 


HI» 如 图 4 


乘法 ， 提 取 谱 域 中 对 


色 点 )， 最 终 ， 将 对 称 中 点 转换 到 空域 ， 进 行 
， 获 得 模型 的 内 列 对 称 平 


中 点 及 平 再 


征 向 


量 的 距离 


称 中 点 (如 


N 


(b) 所 示 。 


(c) 不 同 模型 的 内 下 对 称 平 


面 


图 4 a 


摆 及 模型 顶点 的 法 向 量 


的 法 向 量 为 N， 


ds =Vvev: 


分 布 ,建立 


可 效 识别 


， 所 有 模型 顶点 v 


在 对 称 平面 上 投影 的 中 心 点 为 y，, 则 计算 模型 顶点 到 平面 的 


Ny aD 


d, 0 ,顶点 位 于 对 称 平面 P 的 法 方向 同 侧 ; 当 4d, <0, 
顶点 位 于 对 称 平面 P 的 负 法 向 方向 。 图 


5 (b) 为 基于 内 区 对 称 


KER d, <0, TFH 


j 负 法 向 的 顶点 集 。 


(a) 项 点 与 对 称 平面 关系 (b) 对 称 平 


区 


面 必 <0 的 顶点 分 布 


5 模型 顶点 与 对 称 3 


平面 关系 的 判定 
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在 距离 判断 的 基础 上 ， 引 入 标志 参数 Kk ， 用 来 判 
点 位 于 对 称 平面 法 向 量 方向 ， 
定义 为 无 侧 点 (x =1), 位 于 负 法 向 方向 的 顶点 定义 为 右 侧 


别 模型 顶 


点 (KK 二 一 1); K=0 表 示 顶 点 ， 位 于 对 称 平面 上 。 从 而 获得 具 
有 结构 识别 功能 的 模型 ， 表 示 为 
1 left; 
X(v,,«)|v, esR EC=10 (12) 
-] right; 


通过 顶点 集 d, 的 方向 和 顶点 集 在 对 称 平 
向 一 起 来 确定 的 顶点 的 左右 分 布 。 
根据 模型 表示 X(Cw,x) 以 及 获得 的 模型 内 蕴 对 称 点 对 


MAHE K 的 方 


II,, j Viv; ), 


判别 谱 特征 


点 的 左右 结构 ， 提 取 一 致 性 谱 对 称 特 


征 。 如 图 6(b) 所 示 , 黄色 直线 为 谱 特征 中 左右 结构 的 分 界线 ; 
图 6 Co) 为 依据 统一 结构 所 提取 的 一 致 性 谱 对 称 特征 描述 ， 其 
不 仅 有 效 避 免 了 谱 特征 的 左右 结构 翻转 问题 ， 而 且 提供 了 优化 
的 形状 特征 稀 疏 表示 。 

谱 对 称 特征 提取 算法 如 下 : 


算法 2 求解 3D 模型 的 谱 对 称 特征 描 


5: 


输入 : 3D 空域 自 对 称 点 集合 LL; 。 


输出 : ”3D 模型 的 谱 对 称 特征 。 


a) 依 据 模 型 空域 中 的 自 对 称 点 集合 IILp ， 获 得 谱 域 中 对 称 


点 的 特征 向 量 : 
L= Uola), Q= Qolar): 
b) 计 算 对 称 点 特征 向 量 1; ，Q ;在 谱 域 中 的 欧 氏 距离 


v% e R?” ; 


d(1i,Qj) ， 从 而 获得 谱 域 对 称 中 点 集合 于， 


0) 由 谱 域 中 的 对 称 中 点 于, ， 转 换 为 空域 中 的 对 应 点 ， 并 依 

据 最 小 平面 拟 合 原理 ， 可 获得 空域 中 模型 的 内 列 对 称 平面 ; 
d 由 对 称 平面 以 及 模型 表面 顶点 法 向 量 ， 有 效 识 别 模型 的 

左右 结构 ， 左 向 为 1， 右 向 为 -1， 对 称 面 上 为 0; 建立 具有 结 


构 区 分 的 模型 顶点 表示 ( Xv, K)|v, sR ;Ke(0,-D;); 


e) 由 谱 域 中 的 对 称 中 点 平 , 及 对 应 模型 结构 表示 
X(v ,Kx) ， 获 得 模型 的 一 致 性 谱 对 称 结构 特征 
述 ; 

f)End. 


Re 


(a) 热 核 分 布 。 (b) 谱 对 称 属性 


而 图 7 中 第 一 列 为 2D 谱 特 征 ， 


提取 的 一 致 性 谱 对 称 特 生 
对 称 平面 。 


(c) 谱 对 称 特征 
图 6 不 同 姿态 模型 的 谱 对 称 表示 


BIAR 


EHR, EIAHA 


于 对 称 属性 所 
称 属 性 获得 的 


~ A 
= = 
Li 
图 7 不 同 模型 的 谱 对 称 特征 提取 
2.3 ”基于 一 致 性 点 漂移 的 谱 对 称 形状 配 准 
前 常用 的 配 准 方法 有 和 迭代 最 近 点 算法 E0 (Iterative 


Closest Points, ICP), 薄板 样 条 函数 BA(Thin Plate Spline, TPS), 


多 维 
漂移 算法 (Coherent Point Drift， 
点 集 的 配 准 看 做 
混合 模型 (GMM) 的 形 心 来 表示 ， 


尺度 分 析 B31(Multi-dimension Analysis, 


个 概率 密度 估计 问题 ， 其 中 
Fa ASRS 


CPD), CPD 


MDS) 和 一 致 性 点 


方法 是 把 


通过 最 大 化 似 然 估计 将 GMM 形 
个 点 集 对 齐 


且 使 用 GMM 后 验 概率 获得 两 人 


心 配 准 


到 数据 点 集 


对 应 关系 ， 实 现 非 刚性 模型 配 准 ， 


更 局 的 精 


于 CPD 方法 而 言 , 算法 不 再 是 简单 的 以 两 个 配 准 


， 而 是 通过 算法 建立 其 
的 权 人 1 
， 实 现 非 刚性 模型 


进行 配 ; 


的 非 刚 性 配 准 


oe 其 次 对 
点 列 的 坐标 值 


中 一 个 配 准 点 列 


参数 的 设 定 可 以 进行 点 列 模型 


数据 集 的 高 斯 


的 灵活 性 。 本 文 


基于 模型 的 谱 


对 称 描 
似 性 分 析 。 
算法 3 基于 谱 对 称 的 形状 匹配 


述 ， 进 


输入 : 模型 M4。，M 的 谱 对 称 特 4 


X NxD = [2 Xo al. (每 个 点 Xn ) 为 M, 的 点 集 。 


步 结 合 CPD 方法 ,实现 了 高 效 、 精 准 的 


EFM IK o 


BAB 
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录用 稿 # m, 


Np =1 PLT =Y+GW, 


aur 
o 


Yuxp = [Y1 Y2,- oul (每 个 点 Ym AM, 的 点 iG 
FART Om O) 变换 使 得 yw ERZI Om 0) BX, ZY 1 
一 个 能 一 一 对 应 起 来 。 =D AP DX) -UPX T) + 
输出 : ”对齐 点 列 为 了 ， 对 应 概率 为 乙 。 á 

a) 看 做 混合 高 斯 模型 (GMM) 的 质心 ，x% 看 做 点 的 集合 。 


tr(TTd (PDT), 
b) 求 出 GMM 模型 和 本 ， 就 可 以 得 出 每 个 x 对 应 的 那个 r( DE) 
ym、 Ym BX, PA RATA DN, A ER AI Dian 来 重新 估计 参数 W M o? o 
Plm|x)= POM px |m)/ p(x,) (13) 6) 对 齐 点 集 为 : T=T(Y,W)=Y+GW,。 
| 一 人 通过 P 得 到 对 应 的 概率 P o 
其 中 : P(x | m) = apie XP 267 (Yq WAT (%_ 0)), g)End. 
1 M+l 
p(m) =—» p(x)= 2》,P(m)p(x|m) ， 就 可 以 找 出 一 一 对 应 
M m=1 
we 
c) 估 计 GMM 和 求 出 T(y,,,0) 。 (a) 热 核 分 布 (b) FEF CPD 算法 与 谱 对 称 的 形状 配 准 


本 方法 用 非 刚性 变换 : TI(Y,v)= 了 +v(Y) ， 其 中 


wa - {aa ‘l 


M 
V(Z) = > wiG(z,ym)+W(z) 。 则 上 式 表示 为 | x 
m=] ] 
| j 
T(Y,W)=Y+Gw (14) | | | 
iy) ~、 日 _ 
UP: Gig sop A ， 参 数 是 权 值 9 全。 | i 
dj) 通过 求 GMM BRK HAE Ty, A) 。 通过 
EM 算法 来 估计 参数 。 
M,N 2 
WB: W =0,0? = | I 
m,n=l (c) 蹲 和 跳 匹配 (d) 站 和 蹲 匹配 
初始 化 : 图 8 基于 谱 对 称 的 形状 配 准 


wO0<w<l),fB>0, 2>0, 


: ， 8 为 不 同人 体 姿 态 模型 ， 使 用 CPD 进行 谱 对 称 形状 匹配 
创建 G: ap Ph , 的 实验 结果 。 其 中 ,图 8(a) 为 模型 的 热 核 分 布 ,(b) Ce) (dd 


Sij 
EM Hk, EARI: 为 蹲 和 跳 、 站 和 足 不 同 姿 态 人 体 模 型 的 对 称 性 谱 配 准 。 
(a)E-spet: 计算 了 3 ”实验 结果 与 分 析 
本 方法 的 实验 平台 为 Pentium 3.2 GHz CPU，8G 内 存 的 IBM 
aS PC， 并 使 用 MATLAB R2015b 图 形 可 视 化 软件 和 C++ 等 开发 
= ， 环境 。 
exp 27 |x, -Om + G(m,.)w)| 本 算法 实现 主要 由 三 个 步骤 组 成 : a] 模 型 内 蕴 自 对 称 分 析 ; b) 基 
Sap et Oe, w (270°)? M TARRA W PN AERAR P TAT A ERR: 0) 基于 
= l-w N CPD 的 谱 对 称 形 状 配 准 。 实 验 中 ， 采 用 了 180 种 不 同 姿态 的 非 


刚性 模型 ， 选 取 模 型 数据 库 中 任意 两 个 模型 作为 待 测 数据 ， 对 
模型 进行 对 称 检测 与 谱 对 称 特征 提取 ;， 最 终 , 基于 一 致 性 点 漂 
移 算法 进行 对 称 模型 的 配 准 。 
fi (G+ Ao?d(P1)')W = d(P1) PX -Y 3.1 鲁 棒 性 比较 与 分 析 

本 文 引 入 热 核 (HKS) 和 测 地 线 距离 分 析 模型 自身 的 内 莹 对 


这 一 步 用 旧 的 参数 Ww M o? 来 估计 pj 。 
(b)M-step: 


ChinaXiv 合 作 期 刊 
录用 稿 H m, $: 基于 谱 对 称 的 形状 配 准 方法 
称 性 , 并 基于 谱 域 中 的 特征 距离 有 效 提取 模型 的 内 蕴 对 称 平面 ， (a) 热 核 分 布 (b)ICP 与 局 部 相似 性 配 准 
根据 内 药 对 称 平面 及 模型 顶点 的 法 向 量 分 布 ， 建 立 有 效 识别 模 10 ICP 和 局 部 相似 性 匹配 对 比 结 
型 结构 的 表示 与 谱 对 称 特征 描述 .基于 对 称 特征 描述 的 谱 配 准 


从 上 述 实验 结果 中 可 以 看 出 , ICP 在 配 准时 没 考虑 噪声 的 
影响 ， 若 噪声 干扰 很 大 或 者 是 数据 集 是 线性 的 或 者 是 数据 集中 
有 奇异 点 存在 ， 那 么 用 最 小 平方 和 的 方法 解决 配 准 问题 失效 。 
用 局 部 相似 性 分 析 时 出 现 人 体 模 型 的 有 拘 膊 和 腿 部 等 部 位 相似 性 
较 大 ， 出 现 错误 连接 ， 虽 然 该 方法 能 够 实现 稀疏 化 模型 之 间 的 
分 析 ， 但 是 该 方法 不 能 并 不 能 提供 语义 层次 的 模型 局 部 和 整体 
的 对 应 关系 。 

3.2 ”算法 准确 率 和 效率 比较 

为 了 验证 本 文 方法 的 性 能 ， 对 数据 库 中 的 模型 分 别 用 四 种 方法 
(ICP、TPS、MDS、 本 文 方法 ) 进 行 形状 匹配 统计 。 表 2 为 六 组 
模型 在 1580、1 600、3 200 个 点 进行 匹配 所 需要 的 时 间 。 可 见 
本 文 方法 不 仅 运 算 效 率 明显 优 于 基准 方法 ， 可 实现 更 快速 、 更 
准确 的 形状 匹配 。 
本 方法 利用 内 列 对 称 平面 提取 谱 对 称 特征 ， 有 效 避 免 了 特征 空 
间 存 在 的 左右 结构 翻转 问题 ， 谱 对 称 表 示 不 仅 有 效 保留 了 模型 
的 形状 特征 ， 更 能 够 实现 优化 的 稀疏 表示 ， 提 高 配 准 效率 。1 
表 2 可 见 ， 在 时 间 效 率 上 本 方法 是 相对 突出 的 ， 不 同人 体 姿态 
模型 之 间 的 配 准 率 高 达 90% 以 上 ， 在 变形 大 幅 变 化 的 四 歧 动物 
中 也 具有 85% 以 上 的 配 准 率 (因为 四 肢 动物 本 身 的 结构 性 ), 体 
岗 了 本 方法 在 模型 识别 与 分 类 应 用 中 的 有 效 性 。 同 时 ， 在 与 其 
也 三 种 配 准 方法 比较 中 ， 也 具有 明显 的 优势 。 


优化 了 模型 表示 ,实现 了 精确 、 高 效 的 形状 匹配 .图 9 为 TOSCA 
数据 库 中 部 分 模型 之 间 谱 对 称 的 匹配 结果 。 使 用 本 文 方法 和 
TPS 方法 在 四 组 变形 模型 上 进行 形状 匹配 对 比分 析 , 可 以 看 出 ， 


本 文 的 匹配 结果 明显 优 于 TPS, 对 模型 的 非 刚性 变化 更 加 鲁 棒 。 
图 9 (a) 为 不 同 姿态 模型 基于 热 核 的 分 布 ，Cb ) 为 CPD 对 称 谱 
匹配 结果 ，(c) 为 TRS 对 称 谱 匹 配 结果 。 


l p ICP 是 通过 就 近 点 之 间 的 最 小 距离 约束 模型 的 配 准 ， 若 噪声 干 
(a) 热 核 分 布 (b) 本 文 方法 的 配 准 (CTPS MIB? Bk a Se 
eee ， 扰 非 常 大 或 者 数据 集 是 线性 或 者 数据 集中 有 奇 点 存在 ， 那 么 该 

图 9 本 文 方法 和 TRS 谱 匹 配对 比 结果 a = i z E E P 
方法 解决 配 准 问题 失效 .MDS 是 在 一 定 测度 空间 计算 高 维 数据 
从 图 9 的 结果 可 以 看 出 ，CPD 算法 要 求 待 配 准点 集 (也 ”的 距离 矩阵 ， 然 后 通过 等 距 变 换 将 距离 矩阵 栓 入 到 低 维 欧式 空 


即 高 斯 分 量 的 中 心 ) 作 为 整体 按 给 定 的 几何 变换 关系 进行 一 致 ” 间 中 。 但 其 对 于 非 线 性 的 特征 提取 有 一 定 的 局 限 性 。TPS 是 


| 沪 


Æ 性 移动 ， 因 而 其 能 够 较 好 的 保持 点 集 的 几何 拓扑 结构 ， 即 使 在 ”过 优化 沙 板 样 条 的 光滑 参数 实现 模型 谱 特征 的 非 刚 性 配 准 。 而 
名 不同 姿态 的 干扰 下 仍然 有 较 好 的 配 准 效果 。 本 文 方法 与 其 他 方法 最 大 的 区 别 是 通过 提取 模型 的 结构 性 谱 对 
为 了 进一步 验证 本 文 方法 的 鲁 棒 性 ， 将 本 文 方法 与 ICP MER, 建立 一 个 配 准 点 列 数据 高 斯 混合 模型 ， 用 高 斯 混合 模 

局 部 相似 性 分 析 进 行 了 比较 。 实 验 结果 如 图 10 所 示 。 型 的 质心 代表 原始 点 列 ， 再 与 另外 的 点 列 模型 进行 配 准 。 配 准 

率 不 再 是 通过 求解 距离 和 的 方法 ， 而 是 通过 加 权 GMM 后 验 概 


率 得 到 对 应 关系 。 其 抗 噪声 干扰 能 力 较 强 ， 同 时 提高 了 对 称 匹 
配 的 收敛 速度 ， 配 准 过 程 中 可 靠 性 更 高 。 

图 11 b) Co) 为 使 用 本 文 的 谱 对 称 描述 及 CPD 方法 获得 的 不 
同 模型 之 间 的 形状 匹配 。 不 仅 其 有 效 避 免 了 匹配 过 程 中 的 左右 
翻转 , 对 于 不 同 姿态 模型 的 配 准 及 形状 匹配 具有 较 强 的 稳定 性 。 
图 11 Cc) 展示 了 模型 原始 点 与 对 称 稀 疏 化 的 点 进行 配 准 
的 时 间 对 比 。 
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(a) 四 种 匹配 精度 对 比 (b) 不 同 模型 的 谱 对 称 形 状 匹配 
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(c) 配 准时 间 对 比 图 像 
图 11 不 同方 法 的 精度 对 比 及 不 同 模型 的 谱 对 称 形状 匹配 
4 ”结束语 配 准 。 一 系列 的 实验 验证 本 文 方法 对 于 不 同 姿态 及 不 同 种 类 模 

Za RIA 型 的 识别 具有 较 强 的 稳定 性 和 高 效 性 。 

本 文 针 对 数字 几何 模型 的 数据 量 和 复杂 度 日 meee We 局 特征 进行 对 称 检测 ， 建 
提出 高 效 、 稳 定 的 三 维 模型 形状 对 称 分 析 与 配 准 方法 。 通 过 分 。 立 全 局 的 谱 对 称 特征 描述 比 在 识别 局 部 mens 
析 模 型 热 核 与 几何 上 的 自 对 称 特征 ， ree ete 同时 ”性 方面 还 具有 不 足 之 处 。 
基于 拉 普 拉 斯 映射 的 谱 肉 入 ， 得 出 模型 在 谱 域 中 的 对 称 属性 和 在 未 来 的 工作 中 ， 将 进一步 开展 模型 的 稀疏 化 与 局 部 对 称 
空域 中 的 对 称 平面 ， 有 效 揭示 模型 的 左右 结构 ， 建 立 统一 的 谱 ”性 配 准 研究 ， 建 立 更 高 难度 的 模型 配 准 与 描述 机 制 。 

对 称 特征 表示 ， 进 而 结合 一 致 性 点 漂移 实现 非 刚性 模型 的 形状 
表 2 六 组 对 称 模型 匹配 结果 

模型 原始 点 /个 ” 对称 点 / 配 准 时 间 /s 配 准 率 /% 

di RETE “本文 方 ”本文 方 法 “MDSP IcPl TPSR20 

法 1] 9] 

victoria 1605 582 108. 70 50. 31 92. 96 90. 65 89 93. 12 
〈 站 ) 89 
victoria 
(HE) 1605 175 
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NJ 
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丽 : 对 标的 形状 配 准 束 法 


飞机 1 3200 248 98。23 
飞机 2 3200 305 
鱼 1 1606 463 105。63 
鱼 2 1583 20 
猫 1 3200 908 135。12 
猫 2 3200 194 
狗 1 3200 666 163。16 
狗 2 3200 87 
狗 1 3200 666 113。32 
猫 1 3200 908 
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